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Abstract 

We calculate the spectrum of glueball masses in non-supersymmetric Yang- 
Mills theory in three and four dimensions, based on a conjectured duality 
between supergravity and large N gauge theories. The glueball masses are 
obtained by solving supergravity wave equations in a black hole geometry. 
We find that the mass ratios are in good numerical agreement with the avail- 
able lattice data. We also compute the leading (g^-iV)" 1 corrections to the 
glueball masses, by taking into account stringy corrections to the supergravity 
action and to the black hole metric. We find that the corrections to the masses 
are negative and of order (gy M N)~ 3 / 2 . Thus for a fixed ultraviolet cutoff the 
masses decrease as we decrease the 't Hooft coupling, in accordance with our 
expectation about the continuum limit of the gauge theories. 



*Research Fellow, Miller Institute for Basic Research in Science. 



1 Introduction 



Recently Maldacena formulated a conjecture [1] stating that the large N limit of 
the maximally supersymmetric conformal theories in 3,4 and 6 dimensions are dual to 
superstring/M theory on AdS 4 x S 7 , AdS 5 x S 5 and AdS 7 x S 4 respectively, where AdS d 
is a d-dimensional anti-de Sitter space. More recently Witten proposed [2] that one can 
extend this duality to non-supersymmetric theories such as pure QCD. In this case the 
AdS space is replaced by the Schwarzschild geometry describing a black hole in the AdS 
space. When the curvature of the spacetime is small compared to the string scale and the 
Planck scale, superstring/M theory is well- approximated by supergravity. It was found 
that the supergravity description gives results that are in qualitative agreement with 
expectations for QCD at strong coupling. This includes the area law behavior of Wilson 
loops, the relation between confinement and monopole condensation, the existence of a 
mass gap for glueball states, the behavior of Wilson loops for higher representations, and 
the construction of heavy quark baryonic states [2-7]. 

In this paper, we use the supergravity description of large N gauge theories to compute 
the scalar glueball mass spectrum explicitly for pure QCD 3 and QCD 4 . The glueball 
masses in QCD can be obtained by computing correlation functions of gauge invariant 
local operators or the Wilson loops, and looking for particle poles. According to the 
refinement of Maldacena's conjecture given in [8,9], correlation functions of a certain 
class of local operators (chiral primary operators and their superconformal descendants) 
are related at large N and large (?y M A r to tree level amplitudes of supergravity. The 
correspondence between the chiral operators and the supergravity states has been worked 
out in [8-14]. For example, the operator trF 2 in four dimensions corresponds to the 
dilaton field of supergravity in ten dimensions. Therefore the scalar glueball* J PC = ++ 
in QCD which couples to trF 2 is related to the dilaton propagating in the black hole 
geometry. In particular, its mass is computable by solving the dilaton wave equation [2]. 
In [7], it was shown that the correlation function of Wilson loops is also expressed in 
terms of supergraviton exchange if the distance between the loops becomes larger than 
their sizes, leading again to the supergravity wave equation. 

In this paper we will solve the wave equations numerically to obtain the glueball 
masses. Since this description preserves all the symmetries of QCD, we can identify 
the spin and the other quantum numbers of the glueballs. The mass ratios turn out to 
be in excellent agreement with the available lattice data in the continuum limit. This 

*In the following we will use the notation J p for the glueballs, where J is the glueball spin, and P, 
C refer to the parity and charge conjugation quantum numbers respectively. 
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is surprising since a priori the supergravity computations are to be compared with the 
strong ultraviolet coupling limit of the gauge theory g\ M N 1. 

As we will see, the supergravity computation at gyM-N ^> 1 gives the glueball masses 
in units of the fixed ultraviolet cutoff Aj/y. F° r finite 't Hooft coupling A = g\ M N , the 
glueball mass M would be a function of the form, 

M 2 = /(A)A^. (1.1) 

To take the continuum limit A uv — * oo, we have to simultaneously take A — > so that 
the right-hand side of this equation becomes of the order of the QCD mass scale Aqcd- 
This in particular requires that /(A) decreases as we decrease the 't Hooft coupling A. 

We compute the leading A -1 corrections to the supergravity computation and show 
that this is indeed the case. On the superstring side, the A -1 corrections are due to the 
finite string tension. The leading order string correction to the low-energy supergravity 
action was computed in [15, 16]. This modifies both the background black hole metric and 
the supergravity wave equation in that background. Recently the stringy correction to 
the black hole metric was obtained in [17] by solving the modified supergravity equation. 
We use both this metric and the string corrected wave equation to compute the leading 
A -1 corrections to the ++ glueball masses in QCD 3 . We find: 

(1) The corrections to the masses are negative and of order A~ 3 / 2 : 

/(A) = c + Cl A- 3/2 + • • • , ci<0, (1.2) 

for the ground state and the first 5 excited levels of the ++ glueball. Thus, for a fixed 
ultraviolet cutoff, the masses decrease as we decrease the 't Hooft coupling, in accordance 
with the expectation about the continuum limit of QCD. 

(2) The corrections to the ratios of the glueball masses are relatively small compared 
to the correction to each glueball mass, suggesting that the corrections are somewhat 
universal for all the glueball masses. This may indicate that the good agreement between 
the supergravity computation and the lattice gauge theory results is not a coincidence 
but is due to small A -1 corrections to the mass ratios. 

This paper is organized as follows. 

In section 2, we solve the supergravity wave equations in the AdS§ black hole geometry 
to obtain glueball masses in QCD 3 and compare the results with lattice computations. 

In section 3, we solve the supergravity wave equations in the AdS7 black hole geometry 
to obtain glueball masses in QCD 4 and compare the results with lattice computations. 
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In section 4, we use the string theory corrections to the low-energy supergravity action 
and to the AdS 5 black hole geometry to estimate corrections to the glueball masses in 
QCD 3 . 

We close the paper with a summary and discussions. 

2 Glueballs in Three Dimensions 

The M = 4 superconformal SU(N) gauge theory in four dimensions is realized as 
a low energy effective theory of N coinciding parallel D3 branes. One can construct a 
three-dimensional non-supersymmetric theory [2] by compactifying this theory on R 3 x S 1 
with anti-periodic boundary conditions on the fermions around the compactifying circle 
S 1 . Supersymmetry is broken explicitly by the boundary conditions. As the radius R of 
the circle becomes small, the fermions decouple from the system since there are no zero 
frequency Matsubara modes. The scalar fields in the 4D theory will acquire masses at 
one-loop, since supersymmetry is broken, and these masses become infinite as R — > 0. 
Therefore in the infrared we are left with only the gauge field degrees of freedom and the 
theory should be effectively the same as pure QCD3. 

According to Maldacena [1], the Af = 4 theory in Euclidean R 4 is dual to type IIB 
superstring theory on AdSs x S 5 with the metric 

where l s is the string length related to the superstring tension, g s is the string coupling 
constant and dfl 5 is the line element on S 5 . The x 1,2 ' 3 ' 4 directions in AdS 5 correspond 
to R 4 where the gauge theory lives. The gauge coupling constant g 4 of the 4D theory is 
related to the string coupling constant g s as g\ — g s . In the 't Hooft limit (N — > 00 with 
g\N = g s N fixed), the string coupling constant vanishes g s — > 0. Therefore we can study 
the 4D theory using the first quantized string theory in the AdS space (2.1). Moreover if 
g s N ^> 1, the curvature of the AdS space is small and the string theory is approximated 
by classical supergravity. 

Upon compactification on S 1 and imposing the supersymmetry breaking boundary 
conditions, (2.1) is replaced by the Euclidean black hole geometry [2] 

V - 1 dp 2 + fp 2 - ^) dr> + p 2 f dxj + dfl§ (2.2) 



where r parameterizes the compactifying circle and the a; 1 ' 2,3 direction corresponding to 
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the R 3 where QCD 3 lives. The horizon of this geometry is located at p = b with 

b = w < 2 - 3 > 

Once again, the supergravity approximation is applicable for N — > oo and g s N ^> 1. 

According to [8,9, 12-14], there is a one-to-one correspondence between supergravity 
wave solutions on AdSs x S 5 and chiral primary fields (and their descendants) in the 
Af = 4 superconformal theory in four dimensions. The mass m of a p-form C on the AdS 
space is related to the dimension A of a (4 — p) form operator in the Af = 4 theory by 

m 2 = (A-p)(A + p-4). (2.4) 

The supergravity fields on AdS 5 x S 5 can be classified by decomposing them into spher- 
ical harmonics (the Kaluza-Klein modes) on S 5 . They fall into irreducible representations 
of the 5*0(6) isometry group of S 5 , which is also the i?-symmetry group of the 4D super- 
conformal theory. The spectrum of Kaluza-Klein harmonics of type IIB supergravity on 
AdSs x S 5 was derived in [18, 19]. Among them, there are four Kaluza-Klein modes that 
are SO (6) singlets, coming from the s-wave components on S 5 of bosonic fields. They 
are: 

(1) The graviton g^ v polarized along the R 4 in (2.2). It couples to the dimension 4 
stress-energy tensor of the Af = 4 theory. 

(2) The dilaton and the R-R scalar, which combine into a complex massless scalar field. 
Its real and imaginary parts couple to the dimension 4 scalar operators O4 = tr F 2 and 
O4 = tr F A F of the Af — 4 theory respectively. 

(3) The NS-NS and R-R two-forms, which combine into a complex-valued antisymmetric 
field A^, polarized along the R 4 . Its (AdS mass) 2 = 16 and using (2.4) we see that it 
couples to a dimension 6 two-form operator of the Af = 4 theory. This operator has been 
identified as C 6 = d abc F^F^ 'F$„ [22,23]. 

(4) The s-wave component of the metric and the R-R 4-form polarized along S 5 . 
They combine into a massive scalar with (AdS mass) 2 = 32 and couple to a dimension 8 
scalar operator constructed from the gauge field strength F^ of the Af = 4 theory [24, 23]. 

Only these 50(6) singlet fields are related to glueballs of QCD 3 since 50(6) non-singlets 
are supposed to decouple in the limit R — > 0. 

Let us discuss now how to identify the quantum numbers of the glueballs. The spin 
and the parity of a glueball in three dimensions can be easily found from the transforma- 
tion properties of the corresponding supergravity field. The charge conjugation, C, for 



4 



gluons is defined by A a T^ — > —A^T^ where the T a 's are the hermitian generators of the 
gauge group [25]. In the string theory, charge conjugation corresponds to the worldsheet 
parity transformation changing the orientation of the open string attached to D-branes. 
Therefore, for example, the NS-NS two-form in supergravity is odd under the charge 
conjugation. This is consistent with the fact that it couples to O e , which indeed has 
C= -1. 

From the point of view of QCD 3 , the radius R of the compactifying circle provides 
the ultraviolet cutoff scale. To obtain large N QCD3 in the continuum, one has to take 
g\N — > as R — > so that g^N = g\N/R remains at the intrinsic energy scale of QCD 3 . 
Here g 3 is the dimensionful gauge coupling of QCD 3 . This is the opposite of the limit that 
is required for the supergravity description to be valid. As we mentioned, the supergravity 
description is applicable for g\N 3> 1. Therefore, with the currently available techniques, 
the Maldacena-Witten conjecture can only be used to study large N QCD with a fixed 
ultraviolet cutoff R~ l in the strong ultraviolet coupling regime. The results we find are, 
however, surprisingly close to those of the lattice computation, leading us to suspect that 
(g\N)~ l corrections to the mass ratios are small. In section 4, we will estimate the leading 
(g\N)~ l correction to our computation. 

Consider first the ++ glueball masses. These can be derived from the 2-point function 
of the operator tr F^F^ . In the supergravity description we have to solve the classical 
equation of motion of the massless dilaton, 

% [VgdvQgn = , (2.5) 

on the AdSs black hole background (2.2). In order to find the lowest mass modes we 
assume following [2] that $ is independent of r and has the form $ = f(p)e tkx . Using the 
metric of (2.2) one obtains the following differential equation for /: 

(2.6) 

Since the glueball mass M 2 is equal to — k 2 , the task is to solve this equation as an 
eigenvalue problem for k 2 . In the following we set b = 1, so the masses are computed in 
units of b. If one changes variables to x = p 2 , the equation takes the form, 

d 2 f ( 1 1 1 \d/ k 2 , ft 

+ - + 7 + — TT T-- a , 2 Tn / = °> ( 2 - 7 ) 




dx 2 \x x — 1 x + l/dx 4x(x 2 — 1) 

namely it is an ordinary differential equation with four regular singularities at x — 0, ±1 
and 00. 
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Unlike the equation with three regular singularities (known as the hypergeometric 
equation), analytic solutions are not known for this type of equation. Fortunately there is 
an analytical method to compute its eigenvalues k 2 . It is the exact WKB analysis recently 
developed by mathematicians at RIMS, Kyoto University [20]. To use their approach, we 
note that the differential equation (2.7) can be written as the Schrodinger-type equation 



~ + Q(x)\ g(x) = 0, (2.8) 



where g(x) = ^Jx(x 2 — l)f(x) and 



_ 3a: 4 - 6x 2 - 1 k 2 
Q{X) ~ 4x 2 (x 2 - l) 2 + 4x(x 2 - 1) • { } 



To apply the WKB analysis, one can perturb the equation as 

/ d 2 



dx 2 



+ Q(x) + ( V 2 - l)R(x) g(x) = 0, (2.10) 



by introducing a large parameter 77. With a suitable choice of R(x), the secular equation, 
which determines the values of k 2 so that the equation admits a solution regular at both 
x — 1 and 00, becomes explicitly solvable as a asymptotic power series expansion in rf x . 
Assuming the expansion is Borel summable at r\ = 1, the eigenvalues are approximated 
by the following expression [21] 

k 2 = -6n(n + 1) , (n = 1, 2, 3, ...). (2.11) 

We should note that the differential equation in question is degenerate from the point of 
view of the exact WKB analysis and a mathematical proof of the Borel summability in this 
case has not been given. It is possible that the formula (2.11) receives small corrections. 

Since the analytical expression (2.11) for k 2 is still preliminary and we would like to 
find masses for the other glueball states, we also solved the differential equation (2.6) 
numerically. For large p, the black hole metric (2.2) asymptotically approaches the AdS 
metric, and the behavior of the solution for a p-form for large p takes the form p x , where 
A is determined from the mass m of the supergravity field: 

m 2 = A(A + 4-2p) . (2.12) 

Indeed both (2.6) and (2.12) give the asymptotic forms / ~ l,p -4 , and only the later is 
a normalizable solution [2] . Changing variables to / = ip/p 4 we have: 

(p 2 - p 6 ) r + (3p 5 - 7p) V + (16 + k 2 p 2 ) i> = . (2.13) 
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For large p this equation can be solved by series solution with negative even powers: 

^ = £^ =0 a 2n p- 2 " • (2-14) 
1. The first few coefficients are given 



Since the normalization is arbitrary we can set cto 
by: 



k 2 



fl 2 — 77: 



k 4 



CZ4 = 



a 6 = 



12 
1 

2 + 384 

k 6 

120 + 23040 



For n > 5 the coefficients are given by the recursive relation: 

{n 2 + 4n)a„ = /c 2 a n _ 2 + n 2 a n _ 4 . 



(2.15) 



(2.16) 



Since the black hole geometry is regular at the horizon p = 1, k 2 has to be adjusted so that 
/ is also regular at p = 1 [2]*. This can be done numerically in a simple fashion using 
a "shooting" technique as follows. For a given value of k 2 the equation is numerically 
integrated from some sufficiently large value of p (p ^> k 2 ) by matching f(p) with the 
asymptotic solution set by (2.14) and (2.15). The glueball mass M is related to the 
eigenvalues of k 2 by M 2 = —k 2 in units of b 2 . The results of the numerical work are listed 
in Table I. They agree with the formula (2.11). The 4% discrepancy of the two results 
are either due to some systematic error in the numerical analysis or due to corrections to 
the analytical formula (2.11). 



state 


numerical method 


exact WKB method 


ratio 


0++ 


11.59 


12 


1.03 


0++* 


34.53 


36 


1.04 


Q++** 


68.98 


72 


1.04 


Q-) — |-*** 


114.9 


120 


1.04 


Q-) — |-**** 


172.3 


180 


1.04 


qH — |-***** 


241.2 


252 


1.04 



Table 1: (Mass) 2 of 0++ glueball in QCD 3 obtained by 
solving the supergravity wave equation in the black hole 
geometry (in units of b 2 ) using the two different methods. 



*We thank A. Jevicki and J. P. Nuncs for communication on the boundary condition. 



7 



Since both methods give the same results within a 4% error, we are ready to compare 
them with the lattice gauge theory computations [27]. Since the lattice results are in 
units of string tension, we normalize the supergravity results so that the lightest ++ 
state agrees with the lattice result. The results are listed in Table 2. One should also 
expect a systematic error in addition to the statistical error denoted in Table 2 for the 
lattice computations. 



state 


lattice, N = 3 


lattice, N — > oo 


supergravity 


0++ 


4.329 ± 0.041 


4.065 ±0.055 


4.07 (input) 


0++* 


6.52 ±0.09 


6.18 ±0.13 


7.02 


Q++** 


8.23 ±0.17 


7.99 ±0.22 


9.92 


Q-| — |-*** 






12.80 


Q-| 






15.67 


Q-| 






18.54 


Table 2: 


0++ glueball 


masses in QCD 3 


coupled to 



tr F^F^ . The lattice results are in units of the square 
root of the string tension. The denoted error in the lat- 
tice results is only the statistical one. 

Next we consider the two-form of the supergravity theory. As noted previously, it 
couples to the operator Oq. This operator contains l + ~ and 1 components, which 
correspond to the fields A ri and A^, where i, j = 1,2,3 correspond to the three coordinates 
Xi of R 3 . The remaining components A pT and A pi can be set to zero by an appropriate 
gauge transformation. In the QCD 3 limit R — > 0, the 1 component Aij is reduced to 
a operator in 3D, and thus has a non-zero overlap with the glueball. On the 
other hand, the l + ~ components A Ti couple to an operator which is supposed to decouple 
in the R — > limit. Therefore they do not correspond to glueball states in QCD 3 . 

The s-wave component of the two-form field satisfies the equation [18,26] 

- 16g^g^A M = , (2.17) 

where [ ] denotes antisymmetrization with strength one. As before we look for solutions 
which are independent of r and are of the form A^ = hij(p)e lkx . The pr and the pi com- 
ponents of this equation simply result in a constraint which sets the transverse component 
of A^j to zero. For the remaining pseudoscalar component from A^ the equation reduces 

^The parity P = — 1 is due to the fact that the 2-form is dual to a pseudoscalar. The charge conjugation 
C = — 1 is inferred from the string worldsheet parity. 
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to 

p (p 4 - l) h" + (3 + p 4 ) ti -{k 2 p + 16 p 3 ) h 
in units where 6=1. At large p the solution has the form 

h = p" 4 E~ a 2 „p- 2ra . 



, 



(2.18) 



(2.19) 



Since the normalization is arbitrary we can again set ao = 1. The first few coefficients are 
given by: 

20 

640 + k 4 



a 2 



CL4 



a 6 



960 
3520A; 2 + k 6 



(2.20) 



80640 

We have solved the differential equation and obtained the eigenvalues k 2 by the same 
numerical method described above. The results are shown in Table 3. The supergravity 
results are displayed in the same normalization as the one used in Table 2. 



state 


lattice, N = 3 


lattice, iV — > 00 


supergravity 


0~ 


6.48 ±0.09 


5.91 ±0.25 


6.10 


* 


8.15 ±0.16 


7.63 ±0.37 


9.34 


Q ** 


9.81 ±0.26 


8.96 ±0.65 


12.37 


Q *** 






15.33 


Q **** 






18.26 


Q ***** 






21.16 



Table 3: glueball masses in QCD 3 coupled to Oq. 
The lattice results are in units of square root of the string 
tension. The normalization of the supergravity results is 
the same as in Table 2. 

Since the supergravity method and the lattice gauge theory compute the glueball 
masses in different units, one cannot compare the absolute values of the lowest glueball 
mass obtained using these methods. However it makes sense to compare the lowest glueball 
masses of different quantum numbers. Using Tables 2 and 3, we find that the supergravity 
results are in good agreement with the lattice gauge theory computation [27]: 



^ 0++ ' supergravity 
V M o++ ) lattice 



= 1.50 

= 1.45 ±0.03 . 



(2.21) 
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There are still two more 5*0(6) singlet supergravity fields. One of them is the s-wave 
component of the metric g" and the R-R 4-form A a p jS polarized along S 5 . From (2.4) we 
see that it should couple to a dimension 8 scalar operator 0$- In [23,24], this operator is 
identified as a symmetrized form of F 4 — \{F 2 ) 2 . By using the prescription of Tseytlin 
[28] to symmetrize the group indices, one finds that the operator is even under the charge 
conjugation. This is also seen from the fact that g% is clearly even both spacetime and 
worldsheet parity transformations. Therefore g% has the quantum numbers of the ++ 
glueball. The classical equation of motion of g% is that of a massive scalar with (AdS 
mass) 2 = 32 (in units of b 2 ) [18] on the AdS 5 black hole background (2.2). The mass 
spectrum that we get is given in Table 4. In the g s N — > oo limit the operators O s and trF 2 
are not mixed since they couple to different states in the supergravity theory. However, we 
expect that for finite g s N these operators will mix, thus the full ++ spectrum is expected 
to be given by the interleaving of Tables 2 and 4. For example the ++ ** presumably 
corresponds to the first state in Table 4. 

g% and A a/3l5 



8.85 
12.06 
15.00 
17.98 

Table 4: ++ glueball masses in QCD 3 coupled to 0$, the 
normalization is the same as in Table 2. 

The remaining SO (6) singlet is the graviton g^ v . It couples to the energy-momentum 
tensor T^ u and therefore corresponds to the 2 ++ glueball. It would be interesting to 
compute its mass and compare with the lattice result. 



3 Glueballs in Four Dimensions 

To construct QCD 4 , one starts with the superconformal theory in six dimensions 
realized on N parallel coinciding M5-branes. The compactification of this theory on 
a circle of radius R\ gives a five- dimensional theory whose low-energy effective theory 
is the maximally supersymmetric SU(N) gauge theory with gauge coupling constant 
g\ — R\. To obtain QCD 4 , one compactifies this theory further on another S 1 of radius 
R 2 . The gauge coupling constant g± in 4D is given by g\ = g\jRi = Ri/R 2 - To break 
super symmetry, one imposes the anti-periodic boundary condition on the fermions around 
the second S 1 . 
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According to Maldacena [1], the large N limit of the six- dimensional theory is M theory 
on AdS 7 x S 4 . Upon compactification on S 1 x S 1 and imposing the anti-periodic boundary 
conditions around the second S 1 , we find M theory to be on the black hole geometry [2]. 
To take the large N limit while keeping g\N finite, we have to take Ri <C R.2- In this 
limit, M theory reduces to type II A string theory and the M5 brane wrapping on S 1 of 
radius R\ becomes a D4 brane. The large N limit of QCD 4 then becomes string theory 
on the black hole geometry given by 

" ■ \lr- ■ p ;!2 »r- • // -tlit (3.1) 




P s g!N/An 

with a dilaton ~ p 3 / 4 [29]. The location of the horizon p = b 2 is related to the radius 
R 2 of the compactifying circle as 

As in the case of three dimensions, we will compute the spectrum of glueball masses by 
solving the classical equations of motions of Kaluza-Klein modes of the supergravity the- 
ory. We will consider only singlets of the 50(5) isometry group of S 4 , which corresponds 
to the i?-symmetry group of the six-dimensional theory. 

Consider first the ++ glueball. The non-extremal D4 brane solution has a non con- 
stant dilaton background. As shown in [30] the dilaton is a linear combination of two 
scalars. One of them is massless and couples to the relevant glueball operator. The 
equation of motion for the scalar is given by (2.5) in the background of the metric (3.1). 
Again assuming that the solution is independent of 9 and of the form $ = f(\)e lkx (with 
A 2 = p), one obtains the differential equation in the units where b = 1 as 

(A 7 - A)/"(A) + (10A 6 - 4)/'(A) - A 3 A; 2 /(A) = . (3.3) 

The asymptotic solutions to this equation are / ~ 1, A~ 9 , with the latter corresponding to 
normalizable solutions. In order to solve the equation and find the allowed values of k 2 we 
introduce the function g(X) as /(A) = \~ 9 g(\). This way g(X) has to be asymptotically 
constant for A — > oo, and one can again look for a solution in terms of a negative power 
series in A. The differential equation for g(\) is 

(A 8 - \ 2 )g" + (14A - 8AV - (X 4 k 2 + 5A)g = 0. (3.4) 

The first few coefficients in the power series solution g = J2™=o o-2n^~ 2n are given by (for 
a = 1) 

k 2 

Ct2 — — 

22 



11 



Gt4 



k 4 
1144 
61776 



// 



° 6 " 102960 • (3 ' 5) 
The regularity of / at A = 1, after numerically solving the equation (3.4) as described in 
the previous section, results in the allowed values of k 2 . The first six masses (normalized so 
that the lightest ++ state agrees with the lattice calculation) together with the available 
lattice results [32, 33] are given in Table 5. 



state 


lattice, N = 3 


super gravity 


0++ 


1.61 ±0.15 


1.61 (input) 


0++* 


2.8 


2.38 


Q++** 




3.11 


qH — |-*** 




3.82 


qH — |-**** 




4.52 


qH — 




5.21 



Table 5: Masses of the first few ++ glueballs in QCD 4 , in 
GeV, from supergravity compared to the available lattice 
results. Note that the authors of ref. [33] do not quote 
errors for the ++ * since it is not yet clear whether it is 
a genuine excited state or merely a two glueball bound 
state. 

In order to calculate the masses of the h glueball in four dimensions we will consider 
the 3-form of the eleven dimensional supergravity. In this case, it is more useful to 
use the eleven-dimensional metric 

d\ 2 



ds 2 = 



+ (^-^)^ 2 + A 2 E^ 2 + ^, (3.6) 



(KL - 

U 2 A 4 / 

which reduces to (3.1) upon compactifying x 5 on S 1 and by going to the string frame [2] 
by multiplying the metric by A, setting A 2 = p, and rescaling the other coordinates. The 
s-wave component of the 3-form in the harmonic expansion on S 4 is a singlet of the 5*0(5) 
isometry group [31]. Its mass squared* is 36 in units of b 2 and using (2.4) we see that it 
couples to a dimension 9 operator of the six-dimensional theory. The 3-form obeys the 
following equation of motion: 

\J9 







(3.7) 



*The value of the mass term is fixed by matching the supergravity computation [31] to (2.4) [34]. 
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Choosing a gauge where A pTi and A pi j vanish, where i,j — 1, . . . , 5, and assuming that the 
remaining components are independent of the coordinate r and the only dependence on 
Xi is through e lkx , one finds that there are two independent modes after compactification 
to 4D: 

(1) A three-index tensor A^. This is dual to a massive scalar and can be identified with 
the v glueball of the 4D theory. 

(2) A massive vector A Ti j. This couples to an operator which is supposed to decouple in 
the limit i?2 — *• 0. Therefore it does not correspond to a glueball state in QCD 4 . 

The scalar component of A^ satisfies the differential equation 

(A 7 - A)f (A) + (A 6 + 5)/'(A) - \\e + 36A 2 )/(A) = 0, 



(3.8) 



in the same units as in the equation (3.3) for the dilaton. The normalizable asymptotic 
solution behaves like 1/A 6 , thus we introduce the function g(X) by /(A) = \~ 6 g(\). This 
satisfies 

(A 7 - \)g"(\) - (11A 6 - 17)g'(X) - (72 + A; 2 A 4 )^(A) = 0. (3.9) 
The power series expansion g(\) = S^Lo °2nA~ 2n with a = 1 has the first few coefficients 



a 2 



CZ4 



28 

A: 4 

1792 

129024 + k 6 



fl6 = 193536 ' (3 - 10) 

The single- valuedness of the solution requires g' — 6g at A = 1. With this boundary condi- 
tion, the numerical solution of (3.9) yields the allowed values of /c 2 , and the corresponding 
masses in the units defined above are displayed in Table 6. 



state 


lattice, = 3 


super gravity 


o-+ 


2.19 ±0.32 


1.83 


-+* 




2.67 


Q — h** 




3.42 


q — |-*** 




4.14 


Q 




4.85 


Q 




5.55 



Table 6: Masses of 0~+ glueball in QCD 4 . The lattice 
result is in GeV. 
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Unlike the 3D case, there exists little lattice data on the masses of the excited glueball 
states. We can however compare the ratio of masses of the lowest glueball states CT + and 



and the results are in agreement within the one a error. 

4 Leading String Theory Corrections 

As we mentioned earlier, the supergravity computation is valid in the strong ultraviolet 
coupling limit g s N ^> 1. In order to compare with the lattice computations in the 



so that the scale set by the Yang-Mills coupling constant remains at the intrinsic energy 
scale of QCD. The fact that the glueball masses computed in the supergravity limit are 
in good agreement with the lattice results leads us to suspect that, for this particular 
computation, a' corrections are small. In this section, we test this idea. 

For g s N <C 1, the curvature of the black hole geometry becomes larger than the 
string scale. Therefore stringy corrections (to be precise, the worldsheet sigma-model 
corrections) are expected to become important. The leading stringy corrections to the 
low-energy supergravity action were obtained in [15, 16]. Recently Gubser, Klebanov and 
Tseytlin [17] used the modified action to obtain the leading order string corrections to the 
black hole metric. We use their result to calculate the leading corrections to the glueball 
mass spectrum. We will perform this computation only for the ++ glueballs in QCD3. 
We expect, however, that the conclusions will be similar for the other glueball states. 

According to [17], the leading (in units of the curvature) a' = (AngsN)^ 1 ^ correction 
to the AdS 5 black hole metric (2.2) is 



0++ 





= 1.36 ±0.32 , 



(3.11) 



continuum limit, we have to take g s N — > as we take the ultraviolet cutoff R 1 — > 00 




(4.1) 




(4.2) 
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and 7 is given by 7 = |C(3)o; /3 . With these corrections of the metric, the dilaton is no 
longer constant, instead it is given by 

45 fb 4 6 8 b 12 x 

7 + 27 + V 



There is also a correction to the ten-dimensional dilaton action [15, 16], given by 

l07TGin J \_ I 



(4.4) 

where W is given in terms of the Weyl tensor, and in our background W = 180/p 16 in 
units where 6=1. To the leading order in 7, the dilaton perturbation does not mix with 
the metric perturbation, so we can study the dilaton equation derived from the action 
(4.4) in the fixed metric background (4.1). In subleading order in 7, the term ^e~^W 
would generate a mixing of the dilaton and the graviton. 

We now have all pieces needed in order to obtain the first order correction to the 
dilaton equation. We write $ = $ + f(p)e lkx , with <I> given by (4.3), and expand f(p) 
and k 2 in 7 as 

f(p) = Hp) + iKp) , k 2 = k 2 + ^k 2 . (4.5) 

Here fo(p) obeys the lowest order equation (2.6) and is a numerically given function, and 
/Co are the eigenvalues numerically obtained from the solution of (2.6). The first order 
term of the differential equation obtained from the action (4.4) using the background (4.1) 
and (4.3) is given by 

p -^(^-l) p ^)-klh = (75 - 240p- 8 + 165p- 12 ^ 2/o 



dp \ dp J dp 

+ (75 + 1680p~ 8 - 1815p- 12 )p- lrf/o 



dp 

+ (5k 2 - 120k 2 p' 12 - 405p- 14 )/ (p). (4.6) 

With fo(p) and k% given, one may regard this as an inhomogeneous version of the equation 
(2.6). We solve this equation to determine h(p) and 5k 2 . 

At large p the solution for the first order correction has the form 

h = p~ 4 E~ b 2n p- 2n . (4.7) 

Since (4.6) is inhomogeneous for h(p), one can always add to a solution h(p) a constant 
multiple of the solution /o(p) to the corresponding homogeneous equation (2.6) to obtain 
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another solution. We use this freedom to set b 
given by: 



0. The first few coefficients are then 



b 2 

h 

be 



12 

14400 + 25k 2 k 2 
384 

1344££; 2 + 100800fcg + 35k 2 k* 
23040 



(4.8) 



We can now determine 5k 2 by the same "shooting" method described above for each 
eigenvalue of k% and its corresponding eigenfunction fo(p). It turns out that, for each 
eigenvalue k%, there is is a unique solution with h being regular at p = 1. The first few 
solutions are shown in Table 7. 



state 


(-*2) 


(-5k 2 ) 


Sk 2 /k 2 


0++ 


11.59 


89.75 


7.74 


0++* 


34.53 


365.7 


10.59 


Q++** 


68.98 


809.8 


11.74 


Q-) — |-*** 


114.9 


1397 


12.16 


QH — |-**** 


172.3 


2122 


12.32 


QH — |-***** 


241.2 


2991 


12.40 



Table 7: Leading string correction to the ++ glueball masses in QCD 3 . The 
first column gives the zeroth order supergravity result for the mass squared, 
the second column gives the coefficient of the leading string correction and the 
third column gives their ratio. 

Recalling that the squared mass of each glueball states is given by 

M 2 = - (k 2 + ^A; 2 C(3K 3 + • • •) b 2 , 
we see that the leading stringy corrections to the ++ glueball masses are 



(4.9) 



M 2 ++ 


= 11.59 x (1 + 0.97C(3)a /3 


+ • 


■■)b 2 


Mq+ + * 


= 34.53 x (l + 1.32C(3)a /3 


+ • 


■■)b 2 




= 68.98 x (1 + 1.47C(3)a' 3 


+ • 


■■)b 2 


M 2 


= 114.9 x (1 + 1.52C(3)a' 3 


+ • 


■■)b 2 


M 2 

1K£Q-|--|-**** 


= 172.3 x (1 + 1.54C(3)a' 3 


+ ■ 


■■)b 2 


-/1(/q -)- -)-**** * 


= 241.2 x (1 + 1.55C(3)a /3 


+ • 


■■)b 2 



(4.10) 
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It is important to note that the relation between the radius R of the compactifying circle 
R 4 — > R 3 x S 1 and the location p = b of the horizon also receives an a'-correction. Instead 
of (2.3), we now have [17] 

6= ( 1 _« C( 3 K + ...)_L. (4 . u) 

Therefore, in units of the ultraviolet cutoff A uv = (2/?)" 1 , the glueball masses are ex- 
pressed as 



M 2 ++ 


= 11.59 x (1 


- 2.78C(3)c/ 3 


+ • 


• ")A 2 ^ 


M 2 ++ * 


= 34.53 x (1 


- 2.43C(3)a' 3 


+ • 


• -Wuv 




= 68.98 x (1 


- 2.28C(3)a' 3 


+ • 


■ -)^uv 


M 2 


= 114.9 x (1 


- 2.23C(3)a' 3 


+ • 


■ -)^uv 


M 2 

irjQ-l"-!-**** 


= 172.3 x (1 


- 2.21C(3)a' 3 


+ • 


■ -)^uv 


-A^/q -| — |- * * * * * 


= 241.2 x (1 


- 2.20C(3)a' 3 


+ • 


■ -Wuv 



Thus the glueball masses are indeed modified by the a' = (AngsN)" 1 ^ 2 correction. 
The corrections are negative for all the 6 levels we computed and are of the order a' 3 . 
Therefore the glueball masses decrease as we decrease the 't Hooft coupling A = g s N. As 
we discussed in the introduction of this paper, this is in accordance with our expectation 
about the continuum limit of QCD. 

At the same time, the a' corrections to the ratios of the masses appear to be smaller 
than the corrections to each glueball mass, suggesting that the corrections are somewhat 
universal. This may indicate that the good agreement between the supergravity compu- 
tation and the lattice gauge theory results is not a coincidence but is due to small A -1 
corrections to the mass ratios. Obviously, with the given data, we cannot tell whether 
the stringy corrections for the mass ratios remain small in the continuum limit A — > 0. It 
would be very interesting to see whether this trend continues in the subleading corrections 
in a'. 

5 Summary and Discussion 

In this paper, we computed the glueball masses of large QCD in three and four 
dimensions by solving supergravity wave equations in the AdS black hole geometry. The 
supergravity approximation is valid for large and large A = Qym^ an d therefore the 
results are to be compared with a fixed ultraviolet cutoff in the strong ultraviolet coupling 
regime. 
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We computed the ratios of the masses of the excited glueball states with the mass 
of the lowest state, as well as the ratio of masses of two different lowest glueball states. 
These ratios are in surprisingly good agreement with the available lattice data. We also 
computed the leading A~ 3 / 2 corrections to the glueball masses taking into account stringy 
corrections to the black hole geometry. We found that the corrections to the masses are in 
accordance with our expectation about the continuum limit of QCD. The corrections to 
the ratios of the masses appear to be smaller than the corrections to each glueball mass, 
suggesting that the corrections are somewhat universal. 

The above computations can be generalized to higher spin glueballs. As noted pre- 
viously, the graviton couples to the energy-momentum tensor and solving its equation 
of motion will give the masses of the 2 ++ glueball. In general the higher spin glueballs 
will correspond to operators that couple to massive string excitations. The dimensions 
of these operators are A ~ A 1 / 4 for large A [8]. It would be interesting to see how to 
extrapolate this result to the continuum A — > 0. 

Another interesting issue is the existence of 5*0(6) non-singlet states in supergravity. 
For large A, their masses are of the same order as the SO (6) singlet states we studied in 
this paper. In the continuum limit, Ajjv — > oo and A — > 0, those states should decouple. 
Presumably A -1 corrections make them heavy. 

Maldacena's conjecture reduces the problem of solving large N QCD in three and four 
dimensions to that of controlling the a' corrections to the two-dimensional sigma-model 
with the Ramond-Ramond background. In this paper, we have extracted information 
about glueballs in strongly coupled QCD using the a'-expansion of the sigma-model. It 
would certainly be interesting to understand properties of such a sigma-model better. 
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